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ABSTRACT 

I continue my study of the low-redshift Lya cloud population. Previous work showed how galaxy 
catalogs could be used to attribute relative degrees of isolation to low-redshift Lya clouds found in 
HST/GHRS spectra. This enabled the separation of clouds into two distinct populations corresponding 
to two distinct environments, variously characterized as void/unshocked and non- void/shocked. Void 
clouds have a steep equivalent width distribution [i.e., many smaller absorbers) while non- void clouds 
have a flat distribution. I show that N-body/hydro simulations of Lya clouds are inconsistent with 
observations of the clouds as a function of their environments. Simulations fail to predict the existence 
of significant numbers of detectable void clouds, and incorrectly predict the characteristics of non- void 
clouds. Implicated in this failure is the so-called fluctuating Gunn-Peterson Approximation, FGPA, 
which envisions that Lya absorbers are formed in the large-scale structures of coalescing matter. A 
recent paper (Manning) has modeled the void cloud population as sub-galactic perturbations that have 
expanded in response to reionization. It is notable that success in this modeling was contingent upon 
using the more massive isothermal halo in place of the standard Navarro, Frenk & White, for it was found 
that gravitational restraint on evaporation of baryons is key to producing detectable void absorbers. In 
this paper I extend my modeling of Lya clouds to non-void clouds using the same basic cloud model. 
In the case of voids, clouds are in a quiescent environment, while non- void clouds are thought of as void 
clouds that have accreted to the denser, turbulent intergalactic medium surrounding galaxies, and so 
are subjected to bow shock stripping. Model void clouds are analytically shock stripped, and a column 
density spectrum (CDS) is derived, based on the same halo velocity distribution function as that used 
to explain the void CDS. The non-void CDS produced by shocked sub-galactic clouds are found to 
be capable of producing an excellent fit to the observed non-void CDS without recourse to the FGPA 
mechanism. 

Subject headings: intergalactic medium — quasars:absorption lines - dark matter 



L INTRODUCTION 

N-body simulations now occupy a prominent position in 
the field of astrophysics. Their results promulgate a pic- 
ture of structure formation with which it is hard to argue. 
The reason for this is that, through a variety of means, 
simulations have been designed to be consistent with high- 
redshift Lya forest spectra, and they are integrated into 
the body of what can only be termed the standard model. 
However, because of the way cloud simulations are com- 
pounded, the results at high redshift can hardly be cited 
as proof of the accuracy of their simulations ~ only of their 
self-consistency. 

Rec ent simulations of Lya clouds by iDave fc Tripd 
l)2001|) have extended predictions of models down to a red- 
shift of zero. On the basis of their simulations, they predict 
a column density spectrum (CDS) with a steep cumulative 
spectral slope^, Ssim — —1-15 ± 0.04. Simulations gener- 
ally show that Lya forest absorbers arise from the mildly 
overdense, highly photoionized gas that traces the dark 
matter potentials of large- scale s tructures surrounding 
galaxy concentration s jBi Il993t iMiralda-Escude et al.l 
[1996: Wei nberg et al.Lll997HCroft et alill998HRiediger et all 
1199 8t iDave k. Tripd 120011) In this paper, I wiU generally 



refer to this as the intergalactic medium (hereafter IGM). 
This will be contrasted with the voids, which are not dom- 

^The CDS, as well as the equivalent width distribution function, 
is well-fit by a power-law distribution. Thus the "slope" is the index 
of the power-law. 



inated by galaxies. The mechanism for the formation of 
Lya clouds in these simulations is the fluctuating Gunn- 
Peters on ap p roximation (hereafter F GPA), as explaine d 
in'Bi' (Tm^: "Wch lberg et alJ (IT997I) : ICroft et all lll99l : 
iDave fc Tripp. (£0^" It is based on the assumption that 
gas is substantially stripped or evaporated from all smaller 
halos, forming a roughly homogeneous filament surround- 
ing the fully non-linear coalescing structures. These simu- 
lations have never predicted a significant low-redshift pop- 
ulation of Lya clouds in voids; diffuse gaseous structures 
in voids disperse with the Hubble flow, and today would be 
essentially undetectable ijRiediger et aLllTggsHDave et al.L 
119991 

At low-redshift, galaxy redshift catalogs can be used to 
assess the relative proximity of Lya clouds to galaxies and 
their larger structures. By so doing, it will be possible 
to check on the detailed predictions of the simulations. 
Therefore, the study of low-z absorbers as a function of en- 
vironment may provide a crucial test for the N-body/hydro 
simu lations. 

In lMannind ll2002|) ("hereafter Paper I), low-z Lya ab- 
sorber data frenton et al.ll2000tl was used, , together with 
galaxy redshift catalogs, to assess the relative isolation of 
low-redshift Lya clouds. It was shown how the summed 
scalar tidal field in space could be calculated, so that the 
relative isolation of clouds could be evaluated. I used a 
limiting tidal field Tum to divide the cloud catalog into 
low, and high-tidal field catalogs. As with Paper I, the 
tidal field is given in units of inverse Hubble time squared, 
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Fig. 1. — Panel The cumulative fraction of space in the survey 
IPenton et all 12000 1) with tide T < Tum- Assuming the sightlines 
are randomly selected, this is an estimate of the volume filling frac- 
tion with tide T < Tjim- The three vertical lines correspond to 
the "transition zone" shown in panel B, indicate the possible range 
in value of the void filling fraction. Panel B: The trend in slopes 
(panel a) and intercepts (panel b) of EWDFs defined by catalogs 
with tidal field T < (top of each panel), or T > (bot- 

tom), respectively, for the low-z cloud sample (see §6.2.1 Paper I). A 
strong transition in the slopes in the range —1.3 < logT^^^ < —0.7 
is apparent (two dotted vertical lines); also seen in the intercepts. 
See further discussion in Paper II. 



resulting in a dimensionless parameter of order unity. The 
equivalent width distribution function (hereafter EWDF) 
of these catalogs appear to be well-fit by power-laws. I use 
the model, log dAf{W)/dz = C + S\og (W/63mA), where 
A/'(W) is the number of clouds per unit redshift with rest 
equivalent width > W. Results of the fits can be seen in 
Fig. (see caption). It was a surprise to find that a 
large population of clouds exist in extreme isolation. Void 
clouds have a steep EWDF {Sy ~ —1-6), while non-void 
clouds have a flat slope {Snv ~ —0.5). Here and elsewhere 
I use the subscripts "V" and "NV" to stand for "void" and 
"non-void" . Thus there exist two separate populations of 
clouds. The trends of fitting parameters with variations 
in limiting tide Tum show that the two distinct types of 
clouds are separated by a transition zone (vertical dotted 
lines), whose width is plausibly caused by measurement 
errors, intrinsic scatter (in the TuUy-Fisher relation), pe- 
culiar velocities of clouds and galaxies, and the spatial 
range a cloud travels during its transition from a void- 
type to a non-void-type galaxy. The dichotomy of types 
that is seen in this data, especially in the difference in 
slopes between void (T < Tum), and non-void EWDFs 
(T > Tiira), is dramatic. There exists a similar dichotomy 
- unshocked versus sh ocke d - found in the simulations 
of iRiediger et all l|1998j) and lCen fc Ostrikeil l|199g() . sug- 
gesting that the phenomenology of non-void clouds may 
be explained in terms of shocks. 

The cumulative distribution of tidal field strengths (T < 
Tiim) over the lines of sight fro m which the cloud data 
were taken ijPenton et al.lf2b0Cl(l is shown in Fig. panel 
A. It is used in conjunction with the range of values of 
Tiim in the transition zone (Fig. panel B) to assess the 
volume filling factor of voids. The center of the transition 
is taken to be the tidal field contour at which the shocked 
andunshocked populations meet. At 2 ~ the void filling 



in substantial agreement with the model of lCen fc Ostrike^ 
(1999) (^ 90%). This is the fraction of the universe oc- 
cupied by void-type clouds - hence the volume containing 
unshocked clouds. Void clou ds have a steep slope that 
, according to the analysis of iMannind l)2003() (hereafter 
Paper II), requires clouds to have flat baryon distribu- 
tions {pb oc r~^-^'^), so that the observed absorption sys- 
tems are detected at relatively large impact parameters 
(eg, (rp) ~ 29 kpc for Nhi = 10^^ cm"^). In the shocked 
regions surrounding galaxy concentrations, such a diffuse 
structure would not be possible if significant cloud motions 
were present. Of course the logic behind expecting cloud 
motions is that the gravitational potentials of coalescing 
large-scale structure are driving them. 
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Fig. 2. — The column density spectrum of the non-void clouds 
(lower heavy line for T > 0.3), shown in comparison to the mean 
CDS (light li ne). The dashed line represents the slope of the CDS 
predicted bv lDave &: Tr ipp (2001) for clouds that are expected to 
be found in regions of high tidal field. The non-void CDS has been 
adjusted to represent its line density to the CDS averaged over all 
space. Over non-void space, the non-void CDS has a ~ 7 times 
larger line density. 

Because the N-body/hydro simulations have uniformly 
predicted that, at least at lower redshifts, Lya absorbers 
occupy regions in proximity to galaxy concentrations, their 
predicted CDS spectral slope should agree with the obser- 
vations of the local non-void clouds, since both are refer- 
ring to clouds in the same locations^. Fig. |21is a compari- 
son of predicted and observed spectral slopes. The dashed 
li ne represents th e slope of the cumulative CDS prediction 
of iDave fc Tripd pOOl.1 with an arbitrary normalization. 
Also shown are CDS derived from EWDFs according to 

^Simulations produce a CDS, while observations produce 
EWDFs. The blanketing that occurs in clouds that are not optically 
thin {Njji > 10^^-^; Fig 9 Paper II) results in spectral slopes of CDS 
being generally flatter than those of the corresponding EWDF. 
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the method of §5.3 of Paper II; the heavy hne represents 
the CDS of non-void clouds (T > 0.3), an d the thin hne 
is the mean CDS. The predicted slope of iDave fc TrinrJ 
|2001), Ssim ~ —1.15 ±0.04, does not compare well to 
the observed non-void EWDF slope, Snv — —0.5 ± 0.1 
(Table 3 of Paper I), and agrees even less with the CDS 
slope ~ —0.36. Analysis shows that void clouds are re- 
sponsible for the steepness of the mean EWDF at z ~ 0; 
low EW clouds are quite rare in the turbulent IGM, as 
the cut-off in the non-void CDS in the Fig. [3 shows. But 
simulations have "detected" low-column density clouds in 
the IGM where observations suggest that their detection 
would be very difficult. This suggests that the analysis 
of Doppler parameters in standard N-body/hydro simula- 
tions is incorrect. I address this problem in more detail in 

SI 

The need to explain the void clouds, even while the stan- 
dard model predicts no significant numbers, drives one 
to re-consider the predictions of the "fluctuating Gunn- 
Petersonapproximatioii " (FGPAl ^lBiLll993^l^yeini3erg et all 
Il997t IC^ofteTall ITool IPave fc THpnl l200ir) . But S2 of 
Paper II clearly showed that the predicted Doppler pa- 
rameters for clouds in voids produced by this scenario 
exceed by far the ob served 6- value s of void clouds from 
HST/GHRS spectra llPenton et all |2000). Void clouds 
have Doppler parameters roughly half that of non-void 
clouds (Paper I) . Clearly, this shows that void clouds can- 
not be expanding with the Hubble flow. If they are not 
produced by a fluctuating Gunn-Peterson effect, then how 
else does one explain void clouds except as discrete - i.e., 
as subgalactic structures; remnant survivors from the re- 
ionization epoch? 

In Paper II void Lya clouds were modeled as gas as- 
sociated with sub-galactic halos that has responded to 
the epoch of reionization by evaporating from their ha- 
los, albeit at a gravitationally restrained rat e. A one- 
dime nsional Lagrangian hydro code (,Thoul fc Weinberd . 

was used to track the evolution of the baryons fol- 
lowing reionization at high resolution down to a redshift 
z = 0. The products of these simulations are used to cal- 
culate model CDS. Two d ifferent untruncated halo models 
are used, the lNavarro~al. (1996.. . 1997,) (NEW) halo, and 
an isothermal mass distribution {p oc r~^). 

The NEW halo is the outcome of numerical simulations, 
and would seem to be the first choice for a model. How- 
ever, in the analysis of Paper II, the NEW halo failed 
to restrain the evaporation of baryons sufficiently to pro- 
duce measurable column densities in a void environment 
ai z = 0. However, the isothermal halo proved viable as 
long as the distribution of cloud halo circ ular velocities 
were s t eep; consisten t with that derived bv iKlvpin et alJ 
l|1999() : lKlvph] l|2002D for isolated halos (see §6.2 of Paper 
II). 

According to N-body simulations, an isolated galaxy 
should have a mass distribution in agreement with the 
NEW halo; a steep r^^ density profile for r ^ rmax- How- 
ever, this conjecture is inconsistent with the findings that 
satellite galaxy distributions in groups and clusters follow 
an approximately inverse square num ber density relation 
with radiu s far beyond R max or i?,„v l)Seldner fc Peeblesl 
1197 7: Mc Kav et al.L l200'3) - and have velocity dispersions 
that arc flat to similar distances fZaritskv fc White", '199^ 



IZaritskv et al....l997;.Zabludoff fc Mulchaev..2000;.McKav et 



12002(1 . These resuhs are indicative of an isothermal, rather 
than an NEW mass distribution. 

Recent re-thinking of the McKay et al. results l|Prada fc colleague 
1^03} suggest the above results are affected by interloper 
galaxies which are randomly distributed in space (hence 
in velocity or redshift). However, that small galaxies are 
relatively unclustered, does not imply that they are ran- 
domly distributed: their conjecture is in strong contrast 
with findin gs that satellites hav e a number distribution 
-2-1 llMcKav et all l2002(l around their primaries. 
Prada et al. claim also to have found a method to detect 
interlopers placed into simulated data that can then be ap- 
plied to real data. It is odd that though interlopers are to 
be randomly inserted, their Fig. 3 shows a distinctly non- 
random distribution of interlopers. One wonders whether 
the algorithm for identifying interlopers has been applied 
in a way to pick out things inconsistent with the halo 
model. These doubts prompt the question of whether they 
merely "discovered" what they had presumed at the out- 
set. This issue remains to be fully resolved, but for the 
present, I will consider the original work (McKay et al.) 
to be plausible, and probably correct in its essentials. In 
support of this are two points. First, this assumption is 
at least consistent with the results of Paper II that NEW 
halos were not good cloud models, and that isothermal 
halos are. Second, this particular assumption regarding 
the outer regions of galactic halos is not essential to the 
analysis, but it will be seen that it is consistent with other 
data to be developed in §4.3. 

This same isothermal model is the basis of a successful 
explanation of the void EWDF. Perhaps non-void clouds 
may be explained by a similar approach. 

The goal of the current paper is to explain the nature 
of non-void clouds. I do this in terms of the environ- 
ment in which they are found. I here describe a theory in 
which non-void clouds are produced by discrete clouds. At 
the epoch of their formation, sub-galactic halos must have 
been distributed fairly evenly in the universe - approxi- 
mately half in mildly overdense regions, and the other half 
in mildly underdense regions. As the universe evolved, the 
halos in overdense regions must be carried with the flow 
into the growing zones of shocked gas, as visual presenta- 
tions of structure formation c onvincingly sho w. Thus, ac- 
cording to Birkhoff's theorem l'Birkhofl ll923() . underdense 
regions came to have a locally higher Hubble constant, 
which tends to promote the deepening of the underdensity, 
dispersing their halos and suppressing hierarchical growth. 
Of the halos that are now in non-void space, many may 
have arrived during the primary and secondary infall of the 
coalescing large-scale structure. At later times, whe n voids 
are w ell-established, the "gravitational-repulsion" l|PiranL 
Il997|l of the contents of voids requires that void clouds 
will be ejected at velocities dependent on void sizes and 
the values of the expansion parameter in the void. 

Thus, the above reasoning strongly suggests that the 
great majority of current non-void clouds were, at one time 
or another, standard unshocked void clouds. Hence non- 
void clouds can be characterized with roughly the same 
halo velocity distribution function (HVDE) as derived for 
void clouds (Paper II), though the normalization will be 
different (see below). 

The effect of the accretion of delicate void clouds to 
all he shocked non-void space must be dramatic. When the 
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cloud encounters the dissipated gas in non-void space, the 
tenuously held gas is shock stripped, truncating the cloud, 
resulting in a modification of its EWDF. It is these effects 
that I intend to follow in this paper. The goal is to see if 
void clouds [i.e., as sub-galactic structures), transformed 
by plausible shocks in the intergalactic medium, can be 
used to explain the non-void EWDF at low-redshift. 

The cosmology assumed in this paper is that of Papers 
I and II ~ a standard flat lambda model with h = 0.75. 
The total matter density of voids is either referred to as 
— Vv / Pcrit, or as Oy/rim, where fi™ = 0.3 is as- 
sumed. As in the other papers in this series, I assume 
the cosmic baryon density is ^b/^m = 0.1. Recent re- 
sults from the Wilkinson Micr owave Anisotropy Probe 
(WMAP) l|Spergel fc colleaguest [200,*^ suggest a substan- 
tially larger value ^b/^m = 0.166. The effects of the larger 
value will be noted. The large discovered line density and 
small Doppler parameters of void clouds (Papers I and 
II) strongly suggest that they are to a significant degree 
self-gravitating and discrete. I treat them as such herein. 

I begin the analysis in by studying the effects of 
shocks on model void clouds as a function of their ve- 
locity relative to the IGM. The non-void column density 
spectrum is modeled in fjJI In I discuss some of the 
broader issues touched upon in this paper, and summarize 
my findings in §5. 

2. CONSTRAINING CLOUD VELOCITY 

I model all Lya absorbers as initially taking the form of 
void clouds, and therefore to be consistent with the models 
produced in Paper II. I view the difference in slopes be- 
tween non-void EWDFs and void EWDFs as attributable 
to the relatively extreme conditions of the IGM - a denser, 
hotter, and turbulent environment. I propose that the pri- 
mary factor in the modification of cloud characteristics is 
ram pressure from the relative motions of clouds through 
the IGM. 

To constrain the rate of motion in the IGM to which 
clouds are subjected, I model the effects of shocks pro- 
duced by clouds moving at a range of velocities within 
a medium of density p^v- What ram pressure will ex- 
plain the change in shape and relative normalization of 
the CDS? As noted above, this methodology assumes that 
the transition from void to non-void cloud characteristics 
is currently brought about by the impact of the clouds 
with the dissipated gas near the T = 0.1 contour. 

2.1. Shock Stripping of Void Clouds 

IMurakami fc Ikeuchil ljl994D (hereafter MI94) probed 
the effects of the stripping action of blast waves on mini- 
halos. Mini-halos {e.g.. Roes, 198 a) were once utilize d to 
expla in the Lya forest IjReesL Il986l: IMurakami fc Ikeuchl 
and are closely related to those sub-galactic halos 
studied in the present work. MI94 envisioned two phases 
of stripping - the first being the immediate result of the 
onset of the shock, and the second being the gradual strip- 
ping of the remnant by the relative velocity through the 
medium. MI94 found that the isothermal mini-halos suffer 
significant erosion during the initial shock only where the 
internal pressure of the cloud is less than the ram pressure. 
For the case of a cloud with velocity Vd, encountering a 
dissipated medium with gas of density pg the ram pressure 



on the cloud is, 

Pram=Pgvli. (2) 

MI94 showed that mini-halos could withstand the sec- 
ondary stripping for a sustained amount of time if the es- 
cape velocity from the cloud were in excess of the oblique 
flow within the shock. Some clouds were seen to survive 
in excess of 3 Gyr. 
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Fig. 3. — The density profiles of grown halos processed according 
to the prescription of Eq. Ofor various velocities v^i (as noted in 
the panels). The dotted line represents the hydrogen density in the 
background (critical density) medium. When the cloud halo circular 
velocity Vc = 10.6 kms"-*^, the cloud is totally stripped even at v^i = 
25 kms~^. Larger clouds (higher Vc) survive larger ram pressures. 
In the final panel, cross-section of the Vc = 30 kms~^ cloud when 
Vi-i = 200 kms^^ is only about 1/4 of that when v^i = 100 kms"-*^. 

The initial cloud models are produced according to the 
methods explained in §6 of Paper II. In these simulations, 
200 baryonic bins are used. The flrst strip, which removes 
gas at pressures less than the ram pressure, can be imple- 
mented by the following transformation of density in each 
Lagrangian bin, 

p(z)^p(z)e-(''— (3) 

where the pressure in the cloud in the ith bin is p{i). 
This formulation produces a fairly abrupt truncation of 
the baryon density where the cloud pressure is less than 
the ram pressure. It is assumed that the density of gas in 
the non-void medium (see Eq. ^ is, 

PNV ~ 2 ^ Pent, (4) 

where pcrit — SH^/SttG, a function of redshift. This value 
is in recognition of the fact that this region is undergo- 
ing a mildly non-linear collapse, and therefore the average 
density should be mildly super-critical - ie, ~ 2 pcrit- 

In Fig. 121 1 show the results of subjecting un-stripped 
model clouds (top line in each panel) to ram pressures of a 
gas of the above density, and outfall velocities 25, 50, 100, 
and 200 kms"^ 
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2.2. Bow Shock Stripping 

When the cloud outfall velocity is greater than the sound 
speed of the ambient medium, a bow shock will result. I 
assume the IGM has an average temperature T/vy = 10^ 
K, consistent with an adiabatic sound speed of Cs ^ 17 
kms~^ . This te mperature is consistent with the simula- 
tion of iDave fc'Tr ipp (2001.) for an over-density ^ '^Pcnt- 
In sustaining the first shock, the cloud has been stripped 
down to a volume within which the pressure is greater than 
the ram pressure of the shock. From this point, a quasi- 
stable situation results; a bow shock is set up which stands 
off some distance from the cloud proper. The ambient 
gas penetrates the shock, and is compressed and heated, 
according to the Rankine-Hugoniot jump conditions, pro- 
ducing a shear layer that meets with the cooler, denser 
medium of the cloud. There is a high-pressure zone at the 
head of the cloud due to the presence of the moving cloud 
immediately behind it, causing gas to be diverted around 
it. Gas is then accelerated down the flanks of the cloud by 
the pressure differential. The transverse motion over the 
cooler and denser cloud body may cause Kelvin-Helmholtz 
(K-H) instabilities which could strip away the cloud, layer 
by layer. 

For plane shocks in the Mach range of 3 < < 6 
(~ 50 to 100 kms""'^), the velocity inside a plane shock 
has Mach numbers ranging from 0.49 to ~ 0.46, respec- 
tively. Oblique shocks produce higher Mach values inside 
the shock. This fact, plus the pressure gradient along the 
length of the cloud inside the bow shock means that the 
velocity transverse to the cloud may substantially exceed 
Mach 0.5, perhaps by a factor of 1.5 or so. At the same 
time, the contrast between the average cloud density in 
the un-stripped part of the cloud and the gas between the 
cloud and the bow-shock varies from ~ 50 to 64 to 90, for 
Vci = 50, 75, and 100 kms~^, respect ively. According to 
the 1-D analysis of IVietri et al.l l)1997j) , the growth rate of 
K-H perturbations at 0.7 < < 1.0 is essentially zero 
for adiabatic fluids when the density contrast between the 
cloud and the stripping medium is on order 100. These 
considerations suggest that stripping rates are low for the 
outfall velocities we are considering. 

Anecdotal evidence also suggests that clouds may en- 
dure the stripping of the diffuse IGM gas for long periods 
without significant mass erosion. For instance, the very 
existence of the Local Group cloud compact high veloc- 
ity cloud_CHVC 125+41-207 at a distance 50 ^ d ^ 137 
kpc ijBrtins et al.l l2001j) , with its dense core and a long 
tail, strongly suggests it has survived billions of years of 
stripping. For, since CHVCs are found to have an average 
velocity ~ —1 00 kms~^ relativ e to th e barycenter of the 
Local Group ijBraun fc Burtonl Il999j) . this cloud would 
have to have existed for up to ~ 10 Gyr in order to have 
moved a Local Group diameter of ~ 1 Mpc. Since there is 
no conceivable way of creating this cloud within the Local 
Group in recent times, we must conclude that it has been 
traveling through the IGM at a similar velocity for a time 
on order 10 Gyr, retaining much of its mass. Thus this 
and other CHVCs may be ancient objects that have sur- 
vived billions of years of bow shock stripping; they could 
be the dissipatively stripped and compacted an cestors of 
secondary infall, as explained in the scenario of IMannind 
lll99l . 



While it is reasonable to postulate that smaller clouds 
might not endure this stripping as long as larger clouds, 
there is no way at present, short of detailed hydro/gravity 
simulations, to precisely determine the lifetimes of clouds 
under bow shock stripping. However, we may assume 
they all survive on order 10 Gyr, understanding that this 
may over-estimate the quantity of surviving clouds in non- 
void space at the present time. However, since clouds 
with low halo velocities are stripped much more efficiently 
by shocks, large clouds dominate the cross-section of ab- 
sorbers at the velocities relevant in this problem, so that 
small clouds will have only a small effect on results in any 
case. 
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Fig. 4. — The density profiles resulting from the redistribution 
of baryons to an average density within the zone where the cloud 
pressure is greater than the ram pressure, then smoothed at the edge 
with same exponential used in Eq. |2l Sh own are the distributions 
for the range of halo velocities noted in i|2.c{l in ascending order, 
left to right. These profiles provide the fiducial cloud model for the 
calculation of the non-void CDS. 
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2.3. Turbulent Mixing of the Surviving Cloud 

In the globular cluster formation model of lMannin3 
l)1999j) . the internal baryonic density of self-gravitating 
clouds that are stably responding to supersonic winds {i.e., 
clouds are large enough to survive both stages of stripping) 
must increase with time. This is partly due to the deceler- 
ation of the cloud by ram pressure, but is also the result of 
a cycle of pressure-heating at the front-end of the cloud, 
and radiative cooling, as the denser cloud gas is pushed 
along toward the trailing end of the cloud by the shear 
wind. The K-H instabilities may gradually introduce an 
orderly, toroidal convection pattern in the cloud. This 
physical picture may be very similar for void clouds enter- 
ing the non-void environment. The effect of both cloud de- 
celeration (including a displacement of the baryonic cloud 
from its dark halo), and the fluid waves driven down the 
"fetch" of the cloud's edge, is to disrupt the previous den- 
sity profile within the truncation radius of the shock. In 
the simulations of MI94, the central condensation of the 
baryons within the bow shock is quickly lost. There are 
thus good grounds for assuming that the baryons within 
the shock (excepting the shear layer in between) become 
more evenly distributed; cloud containment switches from 
gravitational to pressure confinement when moving into 
no n- void space. 

For a large cloud, a velocity of 100 kms~^ will trun- 
cate the cloud at about 100 kpc (see Fig. I use these 
values as fiducial in the calculation of the time-scale for 
the transformation from a void cloud to a non-void cloud. 
A useful conversion for the velocity is 100 kms^^ ~ 100 
kpc/Gyr. Thus, it would take ~2 Gyr for the cloud to pass 
completely into the dissipated gas of the filaments. Dur- 
ing, and after the process of shock stripping, the baryonic 
cloud is being decelerated by the ram pressure. I calcu- 
late the time required for the baryonic cloud to be slowed 
enough for it to become displaced from the DM by one 
cloud radius. DM plays an important role in maintaining 
the central condensation of the cloud, and when it is gone, 
pressure gradients can quickly redistribute the gas. The 
deceleration caused by ram pressure is, 

a= . (5) 

■nici 

I am interested in when the displacement is one cloud 
radius; rd — 0.5 ai^. Solving, 

t = {'-A J^. (6) 

\Vcl J V iPNV 

Fiducial values for the density of dissipated gas in the 
IGM are, pNV = '^^bPcrit, and pd ~ lOilbPcrit- Thus, 
since rd/vd = 1 Gyr, t w 3.6 Gyr. In 21 I show that 
the characteristic half-width of the transition zone is ~ 
0.5 Mpc (for 0.1 < T < 0.2). For average cloud baryon 
densities <i 20 ^IbPcrit, the transition time-scale t <^ 5 
Gyr. The cloud thus moves less than the half-width of the 
transition zone by the time it is effectively flattened. 

This general picture is confirmed by Fig. 2 of MI94; by 
the time of Figure 2 d, 1.7 Gyr after the initial shock, 
a shear layer envelops a cloud with a flat density profile, 
bounded by sharp density gradients. The loss of the cen- 
tral condensation of this model cloud was accomplished in 
about one half the time estimated above. 



It is thus expected that the clouds in "undisputed" non- 
void space (ie, logT ^ —0.7; outside the transition zone) 
will approximate a flattened, or random mass distribution 
inside the shear layer of the cloud, so that the average 
density in random sightlines is independent of the cloud 
impact parameter. In the transition zone, one may find 
void clouds, clouds in various stages of the process of being 
transformed into non-void clouds. 

To account for this effect in model clouds, I re-distribute 
the mass within the truncated cloud so that the density is 
uniform, but I apply the same exponential factor used in 
Eq. Oto smooth the edges. In addition, the neutral frac- 
tion of hydrogen in the cloud must be "flattened" as well; 
I substitute the mass- weighted neutral fraction within the 
truncated cloud summed over each bin in the truncated 
cloud. Figure ^ shows the density profiles derived using 
this methodology for various outfall velocities Vd. The 
halo velocities corresponding to the various lines in the 
figure are given by 

Vc = 5.31 X IQOOS" kms~\ 

where n = 16 {i.e., 33.5 kms"^) for the far right-hand 
side line"^, n = 15 for the next at 29.9 kms^^, then 26.6, 
23.7, 21.1 kms^^,.., and so on, for n = 14, 13, 12, etc. 

3. MODELING THE NON-VOID CDS 

The CDS is produced using a method very similar to 
that used to model the void clouds and the HVDF to 
produce the void CDS (§6, Paper II). The difference in 
treatment is that there is an extra variable ~ the cloud 
velocity Vd that results in the stripping of clouds. Each 
model CDS is initially normalized using 0y, the void LF 
normalization, and must be adjusted to approximate the 
observed CDS with a "concentration factor" fmuit, which 
is a function of cloud velocity. 

The first step in modeling the non-void CDS is decid- 
ing which density profile of the shocked clouds works best. 
The upper panel of Fig. |S1 shows normalized CDS for a 
shocked cloud {vd — 100 kms^^) in which the density 
profile is unadjusted, as in Fig. O (short-dashed line), one 
in which it is flattened, as in Fig. 0] (solid line), and an 
average of the two profiles (dot-dashed). This average is 
produced by simply averaging the number of H atoms in 
the respective Lagrangian bins of the unflattened and flat- 
tened clouds. Note that the range of legitimacy of this gas 
profile is expected to be confined to the transition zone 
since clouds are expected to be fully flattened by the time 
they emerge from it. These CDS are shown in relation to 
the observed non-void CDS {Tum > 0.1 (heavier jagged 
line). For reference, the void CDS is represented by the 
long-dashed line. The quality of the fit to the CDS is con- 
sistently good over the range 75 Vd <, 200 kms^^. 

These results suggest that the flattened ram-pressure 
stripped cloud profiles described above produce the best fit 
to the non-void CDS. Thus, results confirm expectations. 
I therefore adopt the fiattened profiles as the preferred 
non-void cloud model, realizing it may not work well in 
the transition zone. 

The next step is to constrain Vd- The lower panel of 
Fig. [31 shows normalized model CDS multiplied by the 

''in the modeling of Paper II, this is the largest halo to survive 
without inside-out collapse, and star formation 
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Fig. 5. — Normalized model column density spectra in relation to 
the observed non-void CDS (heavy solid line) and void CDS (heavy 
long-dashed line). The upper panel presents model CDS for unflat- 
toned halos (short-dash; profiles as in Fig. 3), flattened halos (thin 
solid line; profiles as in Fig. 4), and that resulting from a profile 
made by averaging the two (dot-dashed line). In this figure, the 
cloud velocity is v^-i = 100 kms~^. In the lower panel, model CDS 
are presented for various cloud velocities using fiattened cloud pro- 
files (as in Fig. 4). For 200 > v^-i > 75kms~^, these cloud profiles 
provide a good fit to the slope of the observed non-void CDS. The 
normalization factors fmult Ei.re shown in Fig. 151 



adjustable factor fmuit based on cloud velocities of 50, 100, 
and 150 kms~^, shown in comparison with the observed 
non-void CDS, as derived in Paper I. Note that the high 
quality of the fit to the observations shown when Vd — 100 
or 150 kms"""^, is absent for a velocity as low as 50 kms~^. 

Figure El shows the concentration factors fmuit required 
to adjust model CDS, with their void HVDF normaliza- 
tion, to match the observed non-void CDS for a range 
of cloud velocities. Such a concentration factor is the 
natural result of m odels of hierarchical clustering (e.g., 
iLacev fc Coll. Il99,1). 

It is possible to estimate the value of the concentra- 
tion factor independently from the relative values of the 
normalizations of the HVDFs in void and non-void space. 
According to §7 of Paper II, the void CDS can be best ex- 
plained with the HVDF of the "grown" halo cloud profiles, 
with slope parameter a ~ —1.95, and with normalization 
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Fig. 6. — The relation between the cloud velocity Vd, and the 
multiplicative factor fmult necessary to fit the model CDS (made 
by simulating the stripping process) to that of the observed. The 
significance of the lines are explained in the text. 



(j)y ~ 0.06 (/)* (§6, Paper II), where 0* is the normalization 
of the mean luminosity function (LF). With the i3-band 
LF, I use 0* ~ 0.022 h^ Mpc"^ = 0.0093 h% Mpc"^ (Eq. 
36 of Paper II). Knowing the filling factors, we can write 
the equation. 



01/ /y + (t>*NvfNV = 0* 



(7) 



where 4'*nv normalization of the LF in the non-void 

space (again functionally identical to that of the region 
containing shocked gas), with filling fraction /nv, smd fv 
is given by Eq. [Tl(ie., Jnv = 1 - fv)- Thus, 



Jnv 



6.77 



+3.73 
2.95 



(8) 



where the errors are propagated from the range in fv in 
Eq. n Using the void normalization (jiy for the grown 
halos, the expected concentration factor required to arrive 
at the normalization of the non-void population (Eq. [H]) 
would be 



mult 



113 



-1-62 
-49 J 



(9) 



though this is subject to a few caveats (see below). 

At this point it is convenient to consider how things 
would appear if I had used VLh/Vlm. = 0.16, con sistent with 
the WMAP data ijSpereel fc colleaeuesl l200^ . instead of 
0.10 ([Q. I re-ran the simulations discussed in Paper 
II using this new value, and derived the CDS, which is 
consistent with the observed void CDS. The fit was man- 
aged using the same faint-end slope parameter a = —1.95, 
but the normalization was lower; (py ~ 0.02 0*, ~ 1/3 
of the value derived with the lower baryon density. These 
model clouds were analytically shock-stripped as described 
above, resulting in very similar values of fmuit for a given 
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Vci] apparently the larger cloud neutral fraction, and con- 
sequently larger cross-section at a given column density, 
make up for their lower indicative space density in voids 
when flh is larger. However, because (py is now lower, 
Eqs. |S1 and IHI have different values. The non-void nor- 
malization is shghtly higher {<p%y = 7.02^3 08 4'*) and the 
predicted concentration parameter is about 3 times higher 

if mult — 351_]^g4). 

Due to the transferral of matter from void to non-void 
space, the distributions of halo velocities of void and non- 
void clouds are causally linked. If the slope of the HVDF 
does not change going from the unshocked to the shocked 
environments, and if void clouds are indeed the precursors 
of non-void clouds, then the multiplicative factor needed 
to adjust the CDS of stripped clouds to the observed non- 
void CDS is given by Eq. El 

However, E g. |51 should not be accepted uncritically. 
iKlvpinI 1)200 2*) found that sub-halos have HVDFs consis- 
tent with a slope parameter a ~ —1.72, while isolated ha- 
los have a ~ —1.95 (see §6.2 of Paper II). This is suggestive 
of a flattening of the spectral slope in higher-density areas. 
If the observed non-void clouds have average halo circular 
velocities of Vc — 30 kms~^, then for a = —1.72, there 
would be 2.7 times fewer 30 kms~^ clouds in non-void 
space than if a = —1.95. This, in turn, would require a 
reduction of fmuit by a factor ^ 2.7. That is, Eq. IHlshould, 
by rights, have been fmuit = (t>Nv{vc ^ 'iQ)/(t>v{vc ^ 30), 
for the greater effect of stripping on small clouds results in 
the largest clouds providing the bulk of the cross-section 
of stripped clouds. However, we have no direct informa- 
tion on their number density in the IGM, as we do of C* 
galaxies. Though the higher densities of clouds in non-void 
space gives greater opportunity for hierarchical agglomer- 
ation and a consequent flattening of the spectral slope, 
these clouds are by no means generally sub-halos, for the 
overwhelming majority of non-void clouds in this survey 
are far outside the virial radius of any galaxy (see 

On the other hand, it is more plausible that C* galaxies, 
lying in dense zones of the filamentary structures, could 
have enhanced numbers (relative to the initial, but concen- 
trated halo spectrum) due to hierarchical agglomeration 
of smaller halos - especially of sub-perturbations accreted 
during the primary collapse phase at z ^ 3. I consider 
it plausible, therefore, that the fmuit applicable to clouds 
may be up to ~ 3 times smaller than quoted in Eq. |51 

Fig. El shows the trend of fmuit with Vd neglecting this 
possible factor. Over the range of velocities of relevance to 
this study, the concentration factors for the larger baryon 
density (r2;,/f2„j ^ 0.16) are entirely consistent with those 
of the smaller baryon density. The inset shows a set of 
three vertical lines corresponding to the values of fmuit in 
Eq. [51 and imply cloud velocities in the range 42 <^ Vd ^ 
91 kms^^. The dot-dashed line represents the central 
value of fmuit when J7b/f2m — 0.16, implying a central 
value Vci ~ 130 kms""'^, and a range 100 < Vd ^ 150 
kms^^. 

As noted above, the considerations of hierarchical merg- 
ing in dense zones suggest that fmuit , as calculated by Eq. 
El may be less than indicated by Eq. 7 (i.e., 113 and 350 
kms~^ for fife/rim = 0.10 and 0.16, respectively), per- 
haps by a factor on order 3, although the relative isolation 
of most non-void clouds suggest it is a small fraction of 



this. According to Fig. this would imply a significantly 
smaller velocity. For f^h/rim — 0.10, this would imply a 
velocity Vd 50 kms~^. Recall, however, that the qual- 
ity of the fit of the model non-void CDS to the observed 
declines when Vd <, 75kms~^. For the larger Vlb, cloud 
velocity could be expected to drop from ~ 130 kms~^ to 
^ 80 kms~^, and still give a good fit to the observed CDS 
slope. 

The one thing that may argue against this whole logical 
construct is, of course, that non-void clouds can be under- 
stood in terms of a FGPA. I visit this issue in the next 
section. 

4. DISCUSSION AND SUMMARY 

Further discussion in three broad areas is needed to 
tie this investigation together; the FGPA, the transitionn 
zone, and the line-density of absorbers as a function of T 
in non-void space. 

4.1. FGPA 

It was shown in SO] that the predictions of the FGPA for 
the IGM do not appear to agree with observations. The 
presence of many clouds in void space (Paper I) stands 
in contrast with their relative absence under the FGPA. 
Their absorbers are produced in the slowly varying dissi- 
pated gas of the filamentary structure surrounding galaxy 
concentrations. However, the distribution of clouds propo- 
nents predict is inconsistent with the observed distribution 
of clouds in the same non-void locale (see Fig. 2). These 
indicate a significant problem with the standard model of 
low-z Lya clouds, and hence, by extension, perhaps also 
with that of the high-z Lya forest. 

Although the present work has shown that non-void 
clouds may be explained by the effects of shocks on pre- 
viously unshocked void clouds, this does not mean that 
none of the absorbers have their origin in an FGPA. The 
above analysis suggests that if they do exist, the FGPA 
absorber Doppler parameters should be larger than those 
of sub-galactic halos, which are essentially self-gravitating 
and whose integration path lengths are smaller. Similarly, 
it appears likely that the prediction of very low column 
density clouds in the IGM is due to an incorrect assess- 
ment of the Doppler parameters. For among the initial 
assumptions of the FGPA is that of neglecting turbulent 
effects within the line (which is already broadened by rel- 
ati ve velo cities of order 10 to 40 kms~^), as explained in 
l)1993(l . Required integration path lengths are on order 
1 Mpc at z ~ (Paper II). It seems improbable that in a 
turbulent medium there would not be a broadening of the 
spectral absorption lines over distance such as this. 

What sources for turbulence might there be? If the 
"Birkhoff" effect is propelling centrally condensed clouds 
into the IGM, this represents a major source of kinetic en- 
ergy to drive the turbulence. On the high T-side of the 
non-void universe, energy injection comes with superwinds 
from post-starburst galaxies and the like. If turbulent mo- 
tions of ^100 kms^^ (wc/) exists over scales of 200 kpc 
{2rd)^ then over megaparsec scales, a significant turbu- 
lent dispersion of absorption lines on order ^ 100 kms^^ 
would be plausible. Their Doppler parameters are in fact 
large enough to make the lines appear as a continuum de- 
pression. Since FGPA integration path lengths are weakly 
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correlated with column density, FGPA Doppler parame- 
ters should be rather uniformly ^ 100 kms^^. By con- 
trast, the Doppler parameter histogram of non-void clouds 
finds 6 ~ 60 kms^^ with a dispersion of 15 kms~^ or so 
(Fig. 7 of Paper I; Fig. 155 of Paper II). These hist- 
grams show no sign that there is a higher-6 component 
of absorbers. In fact, the 60 kms~^ broadening seems 
appropriate for a centrally condensed cloud undergoing a 
dissipative interaction with the IGM. Thus, there is no 
sign of FGPA clouds at 6 ^ 100 kms~^, the upper-limit 
for the lPenton et all HoO^data. 

4.2. The Transition Zone 

Much of the modeling of non-void clouds depends on 
the physical picture of a sudden transition from a cool dif- 
fuse void to a warm, reltively dense, and turbulent IGM. 
Let us consider this picture in some detail, motivated by 
the following question: why should non-void space have 
a relatively abrupt end, so that clouds are quickly shock- 
stripped? The fundamental fact about the edge of voids is 
that, on the one hand, diffuse, adiabatically cooled mat- 
ter is balanced by a denser and warmer gas on the other. 
Obviously, the pressure in the latter p^y is much greater 
than that in the former pv- The pressure difference will 
cause an expansion of the IGM into the diffuse void gas. 
Alternatively, the fact that the local expansion parameter 
in voids is greater than the mean, suggests that we may 
consider that the cool, diffuse gas is flowing into the denser 
gas at a velocity sufficient to cause stationary shock front. 
The equation for pressure balance is then. 



PWo^t =PNV -pv, 



(10) 



where i;out = {Hy — Ho)Rv- Solving for Uout in Eq- 10, 



Vout 



f PnvTnv 



UruH V pv 



(11) 



The ambient temperature of non-void space is assumed to 
be 10^ K (see i)2.2|l . while the adiabatically cooled void 
space is - 3000 K (Paper II, §4.3 and Fig. 5). We further 
assume the average gas densities are, 

Py =0.1 fib Peril, 

- -on 

PnV ~ Q Pcriti 

where the former suggests the background (unclustered) 
density of baryons is a tenth of the mean, while the latter 
is Eq. 0] Thus, pnv/pv 67 when = 0.3. These 
values imply that 



Vout 



95 kms 



That is, an outfall velocity of Vout ~ 95 kms"""^ will main- 
tain a pressure discontinuity consistent with what was pre- 
sumed here. This accounts for the strong density gradient 
at the boundary between void and non-void space, and the 
apparent rapidity of the onset of shock stripping on clouds; 
it helps confirm the self-consistency of the physical picture. 



4.3. The Non-Void Cloud Spatial Distribution 

I now consider what might be learned from the trend of 
the cumulative line density of Lya absorbers as a function 
of Tiim- The trend in the log of the intercept C in Fig. 

for the non-void EWDFs (solid line, right-hand side of 
upper sub-panel of panel b) appears to be well-fit by a line. 
I find that the trend is consistent with a cumulative line 
density at W = 63 mA of dM{> Tu^)/dz = IS^T/^;^^^. 
The figure shows that this relation is accurate over the 
range 0.1 ^ T ^ 4. The slope may be taken to give in- 
formation about the radial distribution of clouds about 
isolated galaxies, for the steep dependence of tidal field 
strength with distance from galaxies essentially ensures 
that high tidal field regions are close to a strong con- 
centration of mass. Consider an isolated, centrally con- 
densed body. It follows that the tide varies as ^ (Eq. 
175 of Paper I)'*, so that the cumulative line density is 
dAf{< Riim)/dz oc RJil^^, where Rum is the distance at 
which T — Tiira- This implies a differential line density 
d'^N{R)ldzdR oc i?^^-^^. If the radu of clouds within the 
IGM do not vary systematically with tidal field, then the 
number density of clouds is proportional to the differential 
line density. This exponent, —2.18, is quite close to that 
derived for the distribution of s atellite galaxie s arou nd iso- 
lated parent galaxies found bv iMcKav et alJ l)2002 ): n oc 
R^^-^ - a relation valid in the range 133 <^ R <^ 670 h^g^ 
kpc. For a case in which tides are caused by an C* galaxy 
{vc ~ 161 kms-^ Tullv_& Pierce 2000), with Eqs. 14 and 
16 of Paper I, it can be shown that scalar tidal fields in 
the range 0.1 ^ T ^ 4 (inverse Hubble times squared) 

would be produced at a distance 730 <, R <, 2500 h^^ 
kpc. For reference, the high T-side of the transition zone 
(T ~ 0.2 would occur ~ 2 Mpc from an C* galaxy, making 
the half-width of the transition zone ^ 0.5 Mpc. 

The lower limit of the tidal field range is enticingly 
close to the upper limit of the McKay et al. range. If 
these clouds trace mass, then this apparently quite ex- 
tended, though sparse, "halo" of non-void clouds about 
galaxies may represent an extension of the mass distribu- 
tion entailed by the flat velocity dispersions around gi- 
ant field galaxies th at extend out to Rp ~ 670 h.^^ kpc 
ijMcKav et alil2002|) . If this is so, then there is a chance 
to calibrate how cloud number densities trace matter den- 
sity. This may eventually lead to an estimate of the void 
matter density. 

This analysis strongly suggests that the obs erved CHVC 
popu l ation, thought to be d ark matter-held ijBlitz et al.L 
Il999t [Sternberg et al.L l2002() . and with an average Local 
Grou p barycentric motion of — 100 kms~^ (Braun & BurtonL 
119991) , are representatives of sub-galactic halos that have 
accreted to the IGM long ago. They are, I would suggest, 
the closest examples of the discrete clouds which consti- 
tute the non-void column density spectrum. 

5. SUMMARY 

A self-consistent picture has been built that uses shock- 
stripping in the IGM (non-void space) to transform void 

''As noted in §5.1 of Paper I, galaxies are taken to extend to the 
truncation radius Rt, assumed to be 500 h^^ kpc for an C* galaxy. 
In calculating tidal fields, the galaxy mass is summed to Rt, and the 
tide is calculated as though it were a point mass. Thus the above 
method is consistent even within Rt. 
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clouds into non-void clouds. A connection between the 
relative LF normalizations of the HVDFs in void and non- 
void space was used to calculate a "concentration factor" 
that adjusts the CDS of transformed model void clouds for 
the convergent concentration of clouds into the more com- 
pact non-void space of the IGM. Whether the concentra- 
tion factor can successfully explain the observed non-void 
CDS depends on the cloud velocity and the attendant ram- 
pressure stripping. I have shown that non-void clouds are 
consistent with shock stripped void clouds when this cor- 
rection factor is employed. When cloud velocities Vd ^ 75 
kms^^, the shape of the model CDS derived here is in ex- 
cellent agreement with the observations when the baryons 
in shocked clouds lose their central condensation and be- 
come evenly distributed inside the shear layer. Mergers 
in dense regions of non-void space may increase the num- 
ber density of i2*-galaxies, a change which is accompanied 
by a lowering of the "faint-end" slope parameter a. This 
suggests fewer clouds in non-void space (by up to ^ 1/3), 
lowering fmuit, and suggesting lower outfall velocities. For 
the lower baryon density initially assumed for this and pre- 
vious papers, the implied velocity could be less than 50 
kms~^, but does not produce a CDS that fits well with 
observations. However, for the larger baryon density con- 
sistent with the WMAP results, Vd ~ 100 kms~^ appears 
consistent from the standpoint of the quality of the fit to 
the observed non-void CDS, and the implied concentration 
factor fmuit- 

Indications from the systematic velocity of CHVCs, and 
from the requirements of maintaining a strong pressure 
gradient between void space and the non-void environ- 
ment, appears to require a void "outfall" velocity on order 
Vout ~ 100 kms"^ 
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